Toward an accurate Dark Matter power spectrum beyond BAO scales 
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Although there have been various theoretical studies for nonlinear evolution of Dark Matter in 
Newtonian gravity, the main interest has been the nonlinear shift of Baryon Acoustic Oscillation 
(BAO). In this letter, we propose a new approximated expression for nonlinear Dark Matter power 
spectrum applicable for much beyond BAO scales. In particularly, the proposed expression agrees 
with the result of A'^-body simulation with the accuracy better than 2% up to fc = 1.0 [ft/Mpc] at 
z = 3.0. Even at z = 1.0, the accuracy remains within 5% up to = 0.8 [/i/Mpc]. In doing so, 
we proved the Reg PT proposed by Bernardeau et al. (2011) using an approximation for the kernel 
functions used in the Standard Perturbation Theory, and derived a extended version of the Reg PT. 



One of the key quantity in modern cosmology is the 
matter power spectrum because it contains many im- 
portant informations on the evolution of the universe 
as well as of structure formation. Recent observation 
of the Baryon Acoustic Oscillation (BAO) in the power 
spectrum provides us a new method to precisely restrict 
cosmological parameters [1]. 

Recent progress of the cosmological observation greatly 
motivates various theoretical studies for accurate calcu- 
lation of non-linear matter power spectrum. At present, 
the most successful model is given by the Reg PT 
scheme [4, 5]. Especially, Taruya ct al. (2012) [2] gave 
a theoretical prediction which agrees with iV-body simu- 
lation very well on BAO scales with the accuracy better 
than 1% about up to fc = 0.2 [/i/Mpc] at z = 1. 

We have also proposed a new method called the 
"Wiener Hermite (WH) expansion" which gives an accu- 
rate power spectrum over relatively wide range of wave 
number [3]. There, we established the relation between 
the WH expansion and the standard perturbation the- 
ory (SPT), and showed the equivalence between the WH 
expansion and the F-expansion method [2, 4, 5], and de- 
rived the exponential behavior of the power spectrum in 
each order of the WH expansion (F-expansion) in terms 
of the SPT which has been known in the Rcnormalizcd 
Perturbation Theory (RPT) [6-8]. 

The goal of this letter is to present an accurate mat- 
ter power spectrum much beyond BAO scales. In doing 
so, we, for the first time, give a mathematical proof of 
the Reg PT scheme using approximated kernel functions 
in the SPT. Furthermore, the approximated kernel func- 
tions allow us to extend the Reg PT and as a result we 
derive the predicted power spectrum beyond BAO scales. 



First of all, we review the standard perturbation theory 

1 /2 

(SPT). In the situation / = QJi where / = d\nD/dlna 
is the growth function with D, a being the growth fac- 
tor and the scale factor, it is then known that density 
perturbation and the velocity divergence of Dark Matter 



may be expanded as [9], 
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6{z,k) = ^i?"(5„(k), e{z,k) = -aHY,D"9,,{k), 

n=l n—1 

where H is the Hubble parameter and (5„ and 0„ are time- 
independent. The variable k is the Fourier mode. The 
divergence of Dark Matter 6 is defined by 6* = V ■ v where 
V is the velocity of Dark Matter. Then, the n-th order 
solutions of the perturbation theory arc given by. 
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where p[i_„] = Pi + • • ■ + p„ and F„(pi, . . . , p„) = 
^n([Pi, Pn])- The subscript L means the linearized quan- 
tity. The kernel functions F and G are constructed from 
the mode coupling functions a(ki,k2) = (ki+k2)-k2/A:J 
and /?(ki,k2) = |ki -I- k2p(ki • k2)/2fc^fc| according to 
the recursion relation [9] . 

Various modified perturbation theories partially sum 
the infinite order in the SPT, called "resummation" . 
Since the SPT has exact but formal solutions for any 
order of the perturbation theory using kernel functions 
as shown in Eq. (2), any resummation theory should be 
rewritten in terms of the SPT. However, there is no the- 
ory that is proved a complete relation to the SPT except 
for the WH expansion (F-expansion method). In fact, as 
shown in our previous work [3] , the density perturbation 
S is expanded in the WH expansion method as, 
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x41j([ki,k,])i?W([ki,k,]), (3) 

where iJ*-'-* {i = I,-- - ,r} are stochastic variables and 
satisfy the conditions of Eqs.(27)-(30) in [3]. The coef- 
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ficients of the WH expansion (F-expansion) can be de- 
scribed using the kernel functions in the SPT as follows 
(r > 1): 



and we derive the approximated n-th order solution of 
the perturbation theory from Eq. (2): 
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where n factor is multiplied, because we choose one 
Fourier mode pi from n Fourier modes Pi {i = 1, ■ ■ ■ ,n}. 
Since in the linearized theory the density perturbation 
and velocity is related as 



Fr+2n{[kl,K], [Pl, -Pl, P,i, -Pn])PL{Pi)- 
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The linearized perturbation of energy density and the 
linear power spectrum are expressed by (^^(k) = 

(5j^^(fc)ij(i)(k) and PL(k) = [5')^\k)f, respectively, 
where |k| = k. This expression in Eq. (5) means that 
the density fluctuation with the r-the order of the WH 
expansion is the sum of all the {r + 2n) order of SPT. The 
contribution to the power spectrum from Owjj is given by. 



in real space Eq. (10) becomes 
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The coefficients in the F-expansion are simply related to 
ones in the WH expansion, 

4H = '5J''(fci)---'5i''(fcr)FW(ki,...,k,). (7) 

For the divergence of Dark Matter 9, the same relation 
with 5 is satisfied except for replacing of the kernel func- 
tion F to G. 

Notice that once given the relation between the WH 
expansion and SPT, the problem to solve the non-linear 
evolution of Dark Matter eventually reduces to compute 
the kernel functions in the SPT. Therefore, we next con- 
sider the approximation of the kernel functions. For 
|p[ — >■ 0, the kernel functions are approximated as 



and we finally find 
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6{x) = 5n{x) 5l{x - *L(a;)). 



(13) 



n=l 



Thus, the approximation of Eq. (8) extracts the effect 
of the coordinate transformation of the density pertur- 
bation due to the velocity of Dark Matter. This effect 
would cause the shift of the position of the BAO peak, 
but hardly contribute to the non-linear evolution of Dark 
Matter. 

Third, we give a proof of the Reg PT [2, 4, 5]. Using 
Eq. (8), we can derive the approximated solution of i5^^2n 
in Eq. (5) as follows: 
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We can prove this by the induction in n from the re- 
cursion relation of F and G using a(p[i „],p) — > 1 and 

/3(P[l,ri],P) P[l,n] ■ P/2p^. 

Here, we consider the physical meaning of this approx- 
imation in Eq. (8). For |pj| -^0 {i = 2, • • • the 
kernel function F is approximated as 



where cr^ is the dispersion of the velocity coming from 
the effect of the coordinate transformation in Eq. (13) : 
(T^ = j dpPL{p) /6tt. We multiply the factor n to the first 
term because we choose one Fourier mode from n Fourier 
modes. Here, we must subtract (n — l)(5^^^2nlpi,- - ,Pn->o 
because the first term in Eq. (14) integrates the same 
region defined by pi , • • • , p„ n times in the multiple 
integral in Eq. (5). Then, we find 
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Similarly, we can further show the following expression: following approximation for 
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These expressions of Eqs. (15) (16) are the same with 
Eqs.(51) (52) in Bernardeau et al. (2011) [5]. Thus, we 
find that the Reg PT scheme is equivalent to the approx- 
imation of the kernel functions presented in Eq. (8). 

Finally, in this letter, we consider further approxima- 
tions of Eqs. (15) (16) using Eq. (8). In Taruya et al. 
(2012) [2], the first and second orders of the WH ex- 
pansion (F-expansion) are computed by Eq. (16) and 
Eq. (15), respectively, and the third order of the WH 
expansion is evaluated just as tree level without the Reg 
PT scheme (see Eq. (23)-(26) in Taruya et al. (2012) [2]). 
Here, we consider any orders of the WH expansion using 
Eq. (16). When we meet incomputable terms which in- 
clude the higher order of the SPT than the SPT 2-loop, 
we approximate such terms using Eq. (8). 

For the case of r = 1 and r — 2, P^li from Eq. (6) are 
given by 
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In Eq. (17), F^ei and Freg are the same with ones defined 
by Taruya et al.(2012) (see Eqs. (24) and (25) in the 

(2) 

paper). We have the additional terms for PiJ^h which 
are approximated using Eq. (8) in the condition that a 
Fourier mode in (5^jj(ki,k2) is close to zero: ki — > or 
k? — 0. The simialr calculation for r > 3 leads to the 
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All terms except for exp(— A:^(t^)P33 in P^li{r > 3) are 
the additional terms compared to Taruya et al. (2012). 

Thus, we can compute the coefficients of the WH ex- 
pansion up to any order one needs from Eq. (19). When 

(r) 

we summarize Pyj^^i infinite order, we arrive 

at the following result for the approximated full power 
spectrum: 
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Note that the behavior of our solution at low-fc coincide 
with the resuh in Taruya et al. (2012) and the SPT 2- 
loop solution. In fact, we easily see that when we expand 
the cxponetial factors exp(±A:^(T^), Eq. (21) reduces the 
SPT 2-loop solution. Furthermore, the damping factor 
exp(— A:^o'^/2) in Eqs. (15) and (16) which comes from the 
random motions of Dark Matter particles in Eq. (13) is 
canceled by the correlated motion between two points of 
Dark Matter. For example, it is easily visible for the first 
term in Eq. (21) exp(fc^cr^)P^y. Therefore, our solution 
in Eq. (21) does not have the damping behavior unlike 
the result in Taruya et al. (2012). 

In Fig. 1, we plot the predicted power spectra and TV- 
body simulation results (green symbols) devided by the 
smoothed power spectrum without BAO at some red- 
shifts (z = 3.0, 1.0, and 0.35). The blue doted lines are 

the third order of the WH expansion P-^li + PyvH + ^wh 
corresponding to the results in Taruya et al.(2012). 
Strictly speaking, we have the additional terms even for 

(2) (3) 

-PwH ^^^d -^WH but they hardly contribute to the result. 
The black doted lines are P^li {r = 1, ■ • • , 10} in the 
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FIG. 1. The A''-body and theoretical results are plotted for the cosmological models with the WMAP 5yr [10] (Jim — 0.279, 
Ha = 0.721, fib = 0.046, h = 0.701, Us = 0.96 and as = 0.817 ), at some redshifts {z = 3.0, 1.0 and 0.35). We show the 
ratio of the predicted power spectra to the smoothed reference spectra [11], P/Pnw, (Red solid line: the approximated full 
order in Eq. (21), blue dotted: the WH 3rd order, black dotted: each order from 1st to 10th in WH expansion, and green 
symbols: A'^-body results), and the fractional difference between A-body and analytic predicted results, [PNbody — -P]/fNbody. 
The A-body results used in this letter was presented by [12] ^. We used the program presented by Taruya's homepage to 
compute predicted power spectra ^. 

" These results and their initial conditions at ^i^; = 99 were created by a public A'^-body code GADGET2 and the 2LPT code, 
respectively [13]. These Af-body simulations contains 2,048'^ particles and were computed by combining the results with different box 
sizes 2, 048fe-lMpc and 4, 096/i-iMpc, called Lll-Nll and L12-Aril. 
http://www-utap.phys . s .u-tokyo . ac . jp/~ ataruya/regpt_code .html 



order from left to ritht, and the red solid lines are the 
approximated full power spectrum in Eq. (21). On BAO 
scales our results well agree with the A^-body simulations 
as shown in Taruya et al.(2012). Furthermore, the addi- 
tional terms with the higher order of the WH expansion 
produce the nonlinear matter power spectrum on smaller 
scales. In particular, at z = 3.0 the predicted power spec- 
trum agrees with the TV-body result the accuracy better 
than 2% up to fc = 1.0 [ft./Mpc]. Even z = 1.0, our 
solution behaves like the A^-body result and the accu- 
racy is kept within 5% up to A: = 0.8 [ft,/Mpc]. However, 
at z = 0.35 there appears the large difference between 
the predicted power spectrum and A^-body result due to 
the lack of taking nonlinear effects for the density per- 
turbation. This might be corrected if we consider the 
additional terms having the SPT 3-loop order. 

Finally, we mention a problem of our result. The ap- 
proximation for P^\i [r > 2) includes the linear power 
spectrum P]^ in Eqs. (17) and (19), which leads to the 
BAO behavior in high-fc region for example slightly vis- 
ible at A: = 0.2 - 0.5 [h/Mpc] at z = 1.0. However, this 
behavior is not real, because P^li {r > 2) is derived by 
integrating Pl in Eq. (6), and thus the oscillatory be- 
havior should be smoothed out. In order to get rid of 
this unrealistic behavior, we would have to modify our 
approximate method. We leave this as our future work. 

We can also compute the correlation function from our 
result by truncating the WH expansion at a finite or- 
der one needs as shown in [3]. We would be able to 
apply to the calculation of the various effects including 



the nonlinear evolution of the matter density perturba- 
tions straightforwardly: redshift distortion, bispectrum, 
trispcctrum, and bias effects etc. 
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